Understanding spin transport from the motion in SU(2)xU(l) fields 
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Starting from a continuum constituted by charged tops, we formulate the classical counterpart 
of a previously obtained covariant continuity-like equation for the spin current. Such a formulism 
provides an intuitive picture to elucidate the non-conservation of the spin current and to interpret 
the condition for the emergence of an infinite spin relaxation time. It also facilitates the discussion 
on the spin precession in a one-dimensional quantum wire with Dresselhaus and Rashba spin-orbit 
couplings. Furthermore, we derive the diffusion equations for both the charge and spin densities 
and find that they couple to each other due to the Zitterbewegung. 
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I. INTRODUCTION 

Recently, the coherent manipulation of the spin degree 
of freedom in semiconductor, one of the main subjects 
in spintronics [l|, has attracted much attention for its 
potential applications in future information precessing 
and storage technologies. Understandings on the basic 
properties of the spin transport, spin dynamics and spin 
relaxation in semiconductors are required for the design 
of spin-based devices. In this field, the intrinsic spin- 
orbit coupling (SOC) in semiconductors is considered to 
be an effective route for manipulating the spin degree of 
freedom and generating the spin current. The former as- 
pect is involved in the proposal of Datta-Das spin transis- 
tor 0] , while the latter one leads to the numerous studies 
of the spin Hall effect both theoretically [!, 0, H, 0| and 
experimentally @, S, B [M (HI- However, the SOC also 
brings about the non-conservation of the spin current, 
which makes the definition of the spin current cumber- 
some [12 , [TH [TiJ . A covariant form for the continuity- like 
equation for the spin current was given [l2| in the ter- 
minology of SU(2) gauge potentials. It was shown to 
play an essential role in guaranteeing the consistency of 
a generalized Kubo formula for the linear response to 
non-Abelian fields with different gauge fixings [HI]. As 
the work mentioned above are all based on quantum me- 
chanics, one may ask what is the classical interpretation 
for the non-conservation of spin current and whether this 
classical picture can bring new clues to the understanding 
on the spin transport? 

In this paper, we investigate the motion of electrons in 
the presence of SOCs (regarded as SU(2) fields) as well 
as the U(l) electromagnetic fields. Considering a con- 
tinuum constituted by charged tops, we obtain the clas- 
sical counterpart of the continuity-like equation for the 
spin current. It takes the same form as that proposed 
in the view of the quantum mechanics [121 ] and exhibits 
an intuitive picture for the non-conservation of the spin 
current. This classical picture of the motion of electrons 
is the basis of our discussions on the spin transport in 
present paper. The whole paper is organized as follows. 
In Sec. HH the classical counterpart of the continuity-like 



equation for the spin current is formulated. In Sec. IIII1 
we discuss the precession of the spin orientation in a bal- 
listic quantum wire with both Dresselhaus and Rashba 
SOCs. In Sec. IIVI we derive the equations of motion for 
a single spin, from which we can get the condition for 
the emergence of infinite spin relaxation time. Starting 
with the semiclassical Boltzmann equation, we derive the 
diffusion equations for both the charge and spin densities 
in Sec. [Vj Finally, a brief summary is given in Sec. IVIl 



II. CLASSICAL COUNTERPART OF 
CONTINUITY-LIKE EQUATION 

We start with considering a moving top (classical anal- 
ogy of spin) which rotates at a certain rate. A continuum 
constituted by such kind of tops is completely charac- 
terized by a local velocity field v(r,i), a local particle- 
density field p(r, t) together with a local alignment field 
N(r,t) [16j . Hereafter, a letter in boldface denotes for 
a vector in the conventional spatial space with indices 
i,j = 1,2,3 labelling its components, e.g., r = a^e, 
where e% are the bases and repeated indices are summed 
over; while the one with an arrow over head is a vector 
in the spin space, e.g., N = (N X ,N V ,N Z ). The time 





FIG. 1: (Color on line) The left scheme depicts the time- 
evolution of N at a point r; the right scheme illustrates the 
spatial deviation of N by comparing the fields at two neigh- 
borhood points r and r + Ar for which the parallel displace- 
ment is inevitable. 

evolution of N(r, t) is determined by comparing N at 
different times at the same place, i.e., N(r,t + At) — 
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N(r,t) = Co x N(r, t)At, while its spatial deviation is 
determined by comparing N at different places simulta- 
neously, i.e., N(r + Ar, i) - N(r, t) = ^ X iV(r, i)Air 4 . 
The vector fields a; and fii are natural consequences of 
N being a vector with constant module, then we have 



= J(r,i) x N(r,t), 
(U(r,t) x iV(r,*). 



(1) 



By making use of these two relations together with the 
d d ' • 

density conservation — — h — — = 0, we can easily derive 

at oxi 
a continuity-like equation 



,d_ 
^dt 



LOX 



)*+( 



dxi 



\1a X ) Ji 



(2) 



as long as the natural definitions of the spin density 
a = pN and the spin-current density Ji — pvi N are 
employed. Comparing with the quantum mechanical re- 
sults [12|, one can recognize that cj and Qi correspond to 
the SU(2) gauge potentials rjAo and — r\Ai with r\ — K, re- 
spectively. For a two-dimensional electron gas in narrow 
gap zinc-blende III-V semiconductors, the SU(2) gauge 
potentials have been shown 11211 to be related to the 
Rashba [I?) and Dresselhaus jl8l | SOCs, concretely, 



-> 2m , 

A x = ^-{l3, a, , 



A = -^{a, P, 0), Ao=0. (3) 



In terms of these gauge potentials, the SU(2) "electric" 
and "magnetic" fields can be expressed as 



-d Ai - diAo + r]A x Ai 



Bi = tijkdjAk + -j^ijkAj x Ak, 
which provides a "spin-related" force [l2j 
Ti = Si - a ' + tijkJj ■ &k- 



(4) 



(5) 



d 

In the above, we have employed the notations do = — 

at 

d 

and di = — — for simplicity. 

Equations ([lj and (O are the main relations of this 
section which depicts a classical picture for the motion of 
electrons in both U(l) and SU(2) fields. In the following, 
we will carry on the discussions on the spin transport in 
which these results will be employed. 



one-dimensional quantum wire. For example, a spin- 
polarized current is injected at x — and ballistically 
transported through the wire. We consider the case that 
the Dresselhaus SOC is homogeneous while the Rashba 
SOC can be either homogeneous or inhomogeneous. 

a. We first consider homogeneous a, the precession 
of N can be solved analytically and its three components 
are given by 



^/l + 7 2 cos 9 



N x 
N v 

N z = -^-cos(nx + 9), 



sm(nx 



\/l + 7 2 cos ( 



sin(Ka; + i 



1 + 7 2 ' 
1 + 7 2 ' 



cos 9 



(0) 



where N = (Ng,N%,N§) is the initial spin orien- 
tation, 7 = a/ (3, k — 2m^J a 2 + (3 2 jr\ and tan# = 
(A^q — 77Vq )/-\/l + 7 2 A^q . Figure shows the variation 
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FIG. 2: (Color on line) N x , N y and N z are plotted in unit 
of h/2 as functions of x with a — 1CP 12 eV m, /3 = 5 x 
10" 12 eV m in panels (a), (b) and (c). The initial direction 
No is along z direction in panel (a) and N can be antiparallel 
(parallel) to No as marked by the dark (light) green arrow. In 
panels (b) and (c), the initial direction iVo is along y-axis and 
z-axis directions, respectively. N can only return its initial 
direction. In panel (d), a = j3 = 10~ 12 eV m and N can 
rotate to t/-axis direction at a;=0.25 fim, 0.76 /an (marked by 
orange arrows) with No points to x-axis direction. 



III. PRECESSION OF THE SPIN 
ORIENTATION 

As an immediate application of Eq. fl}, we inves- 
tigate the precession of the spin orientation ./V in a 



of N with respect to a;-axis with a = 1CP 12 eV m and 
(3 = 5 x 10~ 12 eV m. The initial direction N affects the 
precession of N. When No is along z-axis direction, TV 
can either precess to the opposite direction of Nq ( e.g., at 
x=0.35 /im marked by the dark green arrow in Fig.^a)) 
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or go back to its initial direction ( e.g., at x=0.7 /jm 
marked by the light green arrow in Fig.[^a)). However, 
when N is along the x-axis or y-axis direction, TV can 
only return to its initial direction and never achieve the 
opposite direction of N , as illustrated in Fig. [DJb) and 
(c). This can be understood as follows. If the initial 
direction is in the plane perpendicular to the revolution 
axis A x , N will precess in this plane and definitely expe- 
rience the opposite direction of N . Besides, the position 
where N rotates to its initial direction is determined ex- 
clusively by k. It can be seen from Eq. ([6]) that the three 
components of N oscillate with the same frequency k. A 
special situation is a = (3 where N x and N y are formally 
equivalent. Thus N can also precess to the y-axis direc- 
tion (re-axis direction) when Nq is in the x-axis direction 
(y-axis direction), as shown in Fig.^d). 

b. We consider an inhomogeneous case a = ax 
which can be realized by tuning the applied gate volt- 
age on the two-dimensional electron gas. Hence, N 
can be solved analytically in form of series expansion, 

°° b 

namely, N x = - k„V— / +2 , N y 



Kb 



b„ 



-x 



n+l 



n - 



1 



N z = 



71=0 
oo 



b n x n with n a — 2ma/rj and 



-K b N y , 



71 = ' 71 = 

ki, = 2m[3/i]. The coefficients are given by bo = Nq, 
b 2 = — (k Nq - - - 3 ^< 

, , K-lb n K abn-2 

and 6„+2 = — 77 r — -, 

+ (n + 2)(n+l) (n + 2)n 

When a is much smaller than (3, the components of N 

can be approximatively written as 



for n = 2, 3 • • • . 



N x = 



N y = 



N z = 



Kb COS ( 



sm(KbX 



sm( 



K b 

N y 

V-COs(K b X + e'), 

COS &' 

N y 

^-sm(K b x + 6'), 

cos 0' 



XCOs(KbX + 9') 



(7) 



with tan#' = —N§/Nq. We plot these three components 
as functions of x in Fig. [3] with f3 = 5 x 10~ 12 eV m and 
a = 10 -6 eV. If the initial direction is in the y-z plane, 
N will almost precess in this plane for small x since the 
vector potential in the £-axis direction A% is much larger 
than that in the y-axis direction A y . As x increases, 
the amplitude of the oscillation of N x becomes larger 
and N tilts out of the y-z plane (see Fig. [3](a)). If the 
initial direction is along the ir-axis direction, N nearly 
keeps pointing in this direction for small x and begins to 
precess as x increases (see Fig. (3Jb)). 

c. We further consider sinusoid-type inhomogeneity, 
a = bsin(qx). One can solve N analytically with the 
help of series expansion in principle, but insufficient in- 
formation can be obtained because the solution cannot 
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FIG. 3: (Color on line) N x , N v and N z are plotted in unit 
of h/2 as functions of x with f3 — 5 x 10 -12 eV m and a = 
10 -6 eV. The initial direction iVo is along z direction in panel 
(a) and along a;-axis direction in panel (b). 



be expressed either in a closed form or in terms of any 
known functions. Now we calculate it numerically. The 
obtained results show that N behaves like a bi-periodic 
function of x, one period is determined by the gauge po- 
tential A x and the other by q, as shown in Fig. 3] with 
-12 



b = 1CT 12 eV m and q = 10 7 
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FIG. 4: (Color on line) N x , N y and N z are plotted in unit of 
h/2 for a — bsin(qx) with q = 10 7 m _1 . The initial direction 
iVo is along z-axis direction in panel (a) and along a;-axis 
direction in panel (b). 



IV. EQUATION OF MOTION FOR A SINGLE 
SPIN AND ITS APPLICATION 



In this section, we focus on the motion of a single 
charged top n = (n x ,n y ,n z ) with constant module n 
which is the classical analogy of an electron. The equa- 
tions of motion for this top under both U(l) and SU(2) 
fields are obtained as 



57(A) - ViAi) x n{i), 



dnlt) 
dt 

dv ■ — • — » 
m— ^ = Si ■ rt(t) + eE t + e ljk v 3 (B k ■ n{t) + eB k ). (8) 

The first equation can be derived from Eq. |T]) by adopt- 
ing N(r,t) = n(t)5(r — r(t)) with f(t) being the trajec- 
tory of a single top. The second equation is due to the 
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fact that its translational motion is governed by both 
the Lorentz force and the "spin-related" force given in 
Eq. ((5} . Clearly, equation @ gives rise to some immedi- 
ate consequences: 

Case 1 The first equation for the time rate of n clearly 
manifests that the vector n does not precess when it is 
parallel to _4o — ViAi, or more specially, Aq — ViAi — 
0. In such cases, the spin orientation keeps unchanged, 
which results in an infinite spin relaxation time. A typical 
example is that the infinite spin relaxation time occurs 
in the ±[1, ±1, 0] direction when Aq = and a — ±/3, as 
discussed in Ref. [l9[ . This is analogous to the case in the 
classical electrodynamics where an electron moving in the 
uniform orthogonal electromagnetic fields with certain 
velocity, v = ExB/B 2 , does not feel the Lorentz force. 

Case 2 In virtue of the coupling between n and the 
SU(2) field in the second equation of ([5]), the time rate of 
ft leads to the time-dependent effective fields even when 
the SU(2) fields are time- independent. 

As a concrete example, we consider a spin Hall system 
with constant Rashba and Dresselhaus SOCs. Thus Aq = 
0, £i = 0, Bi = ^e-ijkAj x Ak and we first neglect the 
electric field. The second equation of Eq. ([5]) reduces to 



m ~dT = ^ ' ( Vj ^ X 

which gives rise to a relation between ft and Vi 

1 



«i(t) 



-Ai-n(t) + C h 



(9) 



(10) 



where Cj = voi + -^Ai ■ uq are determined by the initial 
values i>i(0) = vq{ and n(0) = (ng, 77q, ng). Consequently, 
one only needs to solve one equation, 



^ = 7](—Ai ■ ft - Ci)(Ai x ft), 
at to 

which can be explicitly written as 



(11) 



2m 



dn x 

dt r\ 
dn v 2m 

dt T] 
dn z 2m 



n 



dt 



V 



2m 

1 

rj 



ad - (3C 2 
Pd - aC 2 - 
ad - pC 2 - 
Ci/3 — C 2 a 



AaPn x + 2{a 2 +P 2 )n y 



4a{3n y + 2{a 2 +(3 2 )n x 



4a(3n x 



Aafiny 



(12) 



For a = (3 which is called ReD field in Ref. [l9J], we can 
solve these equations analytically: 



\[2 cos ip 

n" 



1 



y2 cos <p 

n~ 



^—sm{u)t + (p) + -{nl+n%), 



cos ip 



— COs(Ldt + if), 



(13) 



where tan ip 



ng)/(-\/2n§) is determined by the 



initial conditions. It is clear that the tip of 77 experiences 
a cyclotron rotation with frequency w = 2\[2ma{ya x — 
v 0y)/v- The instantaneous velocity solved from Eq. (fTU|) 
is just its initial value v x — vq x and v y — va y , i.e., the 
electron undergoes a motion with uniform velocity. This 
is due to that the time-dependent parts of n x and n v only 
differ from each other by a minus sign. Specially, when 
v x = v y , the spin vector n does not precess since u> = 0, 
which recovers the result in Ref [l9T |. 

When an external electric field E — (E x ,E y ) is ap- 
plied, which mimics to the usual spin Hall effect in cur- 
rent literature, we obtain Vi(t) = —^Ai-n(t) + Ci + ^Eit. 
The equation of motion for 77 is almost the same as 
Eq. ([TTJ) if C t is replaced by C t {t) =C t + ^E t t. For the 
electric field is sufficiently weak, the perturbation theory 
is applicable and we can expand 77 in power series of the 
electric field, 77 = fv a ' + + • • • . The Oth-order results 
take the form of Eq. (TT3|) , while the first order corrections 
are 



77*« = 

y(i) - 



n 



A 1 ) - 



-\{E) t 2 cos(u>t + ip), 
X(E) t 2 cos(ujt + (p), 
-V2X(E) t 2 sm(ujt + (p), 



(14) 



with A(£') = ean,Q(E x — E y )/ (rjcostp). The form of 
the entire motion for 77 is similar to the case without 
external electric fields, but with a time-dependent am- 
plitude d(t) = [(nf ) ) 2 /(2cos 2 (p) + A 2 ^] 1 / 2 and a phase 
(p{t) = Lp + tan _1 (v / 2Ai 2 cos^/ng). 



V. DIFFUSION EQUATIONS AND 
ZITTERBEWEGUNG EFFECTS 

In the presence of U(l) and SU(2) fields 7777^ = 
Pi(v,t) — eAi — r\Ai ■ N, we can write out the charge 
current epvi and spin current pViN as follows 

eji = — [Pi - eAi) At ■ cr, 

m to 



Ji = —{Pi- eAi) ~ —N(Ai ■ a) 
m 777 



(15) 



The first terms on the right hand sides of Eq. (fTS")) cor- 
responds to the conventional non-relativistic currents. It 
is worthwhile to pay attention to the second terms which 
are related to the Zitterbewegung. It was ever believed 
that the Zitterbewegung (20l | as a relativistic effect can- 
not be observed directly due to the high frequency (of 
order 10 20 Hz) and short length scale (of order 1pm) [2l| . 
Recently, J. Schlicmann et al. suggested to detect the Zit- 
terbewegung in III-V semiconductor quantum wells [22| 
where the energy and length scale become available in 
experiments. We will show that the existence of Zitter- 
bewegung-velsded phenomena in the spin transport makes 
the diffusion equations for the spin and charge densities 
couple to each other. 
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In comparison to the expression of the conventional 
non-relativistic current density, there is an extra term 
j'i = x &\i m the expression derived from the Dirac 
equation in the non-relativistic limit. This term 7' can 
be regarded as a result of the Zitterbewegung [2l|, [23| . As 
long as we distinguish between the spin space and the 
conventional spatial space, we have 



7' - —A- 
m 



— [Vp x Nl 

m 



(16) 



In the above derivation, we have adopted Eq. |T]) and the 
fact that the matrix (A)y = A\ is traceless. Hence we 
can consider ^-Ai ■ a arising from the Zitterbewegung for 
either uniform p or Vp cx N. We will see in the following 
that ^Ai ■ a enters into the diffusion equations for both 
charge and spin densities which couple to each other. 

The coupled charge-spin diffusion equations in a spin 
Hall system have been investigated by several groups in 
the view of quantum kinetic theory [24j |. Starting from 
the Boltzmann equation as well as the above results, we 
can derive the semiclassical diffusion equations for both 
charge and spin densities. The introduced distribution 
function /(r, p, t) obeys the Boltzmann equation 



df , df df _df 

n. 1 ' * o — <~\< coll* 

Ot OXi dpi ot 



(17) 



For charged top, Fi should contain both the Lorentz force 
caused by the U(l) fields Ei and Bi as well as the "spin- 
related" force caused by the SU(2) fields £j and Bj. The 
right-hand-side of Eq. (jTTJ) represents the collision term 
contributed by the elastic scattering between electrons. 
For the system not driven far away from the equilibrium, 
the relaxation-time approximation is applicable and the 
distribution function / can be decomposed into the equi- 
librium distribution fo and a small deviation /', in which 
/o is independent of the direction of p. For the spin Hall 
system with constant Rashba and Dresselhaus couplings, 
the magnetic field and the SU(2) "electric" field vanish 
and eijkVjBk ■ N^- is of the 2nd order and can be ne- 
glected. Hence equation (|17l) is written as 



dt mm dxi 



f-fo 



(18) 



~ d p - 

where V, = V, + eEi— with e = and g = fN is the 

oe 2m 
distribution function for the spin density. 



Integrating Eq. (fl8|) over p, we have the diffusion equa- 
tions for the charge density p = e J fdp. In the calcula- 
tion, we encounter J( £i /)dp = f(—f')dp. For a steady 
state, /' can be solved from Eq. (fl8)) 



-r^V 3 f-^A 3 -N^L) 

m m oxj 



(19) 



Neglecting higher orders of /', we have 



/ 



dp (^/') ~ / ,/p/„ = -DY,,>. (20) 



Here D = st is the diffusion coefficient which turns to be 
v^t /2m in the quantum kinetic equation with Vjr being 
the Fermi velocity. As a result, we have 



dp ~ 2 77 - da p-p 

— - DV p Ai ■ — = . 

ot m OXi t 

and the charge current is given by 



Ji 



dp Vi f = -D\7 lP - —Ai 



(21) 



(22) 



We derive the diffusion equation for the spin density in 
an analogous procedure, 

d ~ 
(7^ - riAa y.)a- D(Wi + nAi x )(V< + rjAi x )a 

-^{A i .N)N^ = -^ (23) 
m dxi t 



with the spin current being 



Ji = ~D(Vi 



r/AiX)a -(Ai ■ N)a. 

m 



(24) 



Eq. (|2ip and Eq. (|23|) are classical counterparts of the 
coupled spin-charge diffusion equations [24|. They are 
valid for a large time scale since our startpoint is the 
semiclassical Boltzmann equation which implies the as- 
sumption of the energy conservation. For a short time 
scale St, Ss is large and the quantum kinetic theorem 
should be employed. 

The last terms on the left hand sides of both Eq. (|2"Tj) 
and Eq. (|23|) are the contributions from the Zitterbe- 
wegung which make these diffusion equations coupled 
to each other. To illustrate the effect of the Zitter- 
bewegung-vel&ted contributions, we consider the one- 
dimensional quantum wire with constant Rashba and 
Dresselhaus SOC. The charge density exponentially de- 
cays in the wire, namely, p(x) cx e~~ x l Ls where the de- 
cay length is given by L s = (^ 2 +4D/t + £,/(2D))- 1 
with f = 2(a 2 - P 2 ){aNy + f3N$)/ri(a 2 + /3 2 ), shorter 
than that without the Zitterbewegung-rel&ted contribu- 
tions L s = (V-Dt) -1 . 



VI. SUMMARY 

Depicting moving electrons as a continuum constituted 
by charged tops, we investigated the motion of elec- 
trons in both U(l) and SU(2) fields and discussed the 
spin transport. We derived the classical counterpart of 
the continuity-like equations for the spin current, which 
takes the same form as that ever proposed quantum- 
mechanically in a previous paper [12j. This provided us 
an intuitive picture for elucidating the non-conservation 
of the spin current in the presence of the SOC. We dis- 
cussed the precession of the spin orientation in a ballis- 
tic one-dimensional quantum wire with both Rashba and 
Dresselhaus SOCs and found that the initial direction can 
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greatly affect the spin precession. We also formulated the 
equations of motion for a single spin in the presence of 
both U(l) and SU(2) fields. As a direct consequence, we 
obtained the condition for the emergence of infinite spin 
relaxation time. A special situation of our conclusion 
recovers the previous result discussed by other authors. 
Furthermore, we derived the semiclassical diffusion equa- 
tions for both the charge and spin densities. We found 
that the Zitterbewegung causes these equations coupled 



each other and makes the decay length of the charge den- 
sity much shorter in a one-dimensional quantum wire. 
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